Magnetocaloric effect in spin- 1/2 XX chains with three-spin interactions 
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We consider the exactly solvable spin-1/2 XX chain with the three-spin interactions of the XZX+ 
YZY and XZY — YZX types in an external (transverse) magnetic field. We calculate the entropy 
and examine the magnetocaloric effect for the quantum spin system. We discuss a relation between 
the cooling/heating efficiency and the ground-state phase diagram of the quantum spin model. We 
also compare ability to cool/heat in the vicinity of the quantum critical and triple points. Moreover, 
we examine the magnetocaloric effect for the spin-1/2 XX chain with three-spin interactions in a 
random (Lorentzian) transverse magnetic field. 
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I. INTRODUCTION 

In general, the magnetocaloric effect (MCE) refers to 
any change of the temperature of the magnetic material 
under the variation of external magnetic field. The re- 
vival of interest toward the various aspects of the physics 
of MCE which has been observed recently is mainly con- 
nected with the potential room-temperature cooling ap- 
plications (see Rcfs. 1,2 for recent reviews). Another 
important application of the MCE is the possibility to 
map out the H-T phase diagram by detecting the mag- 
netocaloric anomalies at a magnetic phase transition at 
high (pulsed) fields. For some materials there is no 
alternative way to do that. Since the first success- 
ful experiment of adiabatic demagnetization 3 , MCE is 
the standard technique of achieving the extremely low 
temperatures . 

Another important issue of the MCE is its intimate 
relation with the quantum critical points (QCPs) 5 . The 
MCE can be quantified by the adiabatic cooling rate 
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where Ch is the heat capacity at the constant magnetic 
field, and M is the magnetization. The dependence of the 
cooling rate on the magnetic field is an important charac- 
teristic of a specific magnetic material. The cooling rate 
TTh is related to the so-called generalized Griineisen ra- 
tio, 
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the important quantity characterizing the QCP. It is 
known that the generalized Griineisen ratio changes 
its sign when the parameter r governing the zero- 



temperature quantum phase transitions crosses its criti- 
cal value r c , i.e., in the QCP 6 ' 7 . In the case of MCE r 
in Eq. (1.2) is the external magnetic field H and QCP 
corresponds to the critical value H c at which the sys- 
tem undergoes the transition between different magnetic 
structures at zero temperature 8-10 . As the sign of the 
cooling rate depends on the way magnetic field affects 
the entropy at isothermal conditions, the system can un- 
dergo adiabatic cooling as well as adiabatic heating under 
the increasing (or under the decreasing) of the external 
magnetic field magnitude. Thus, the magnetic materials 
with complicated structure of zero-temperature (ground- 
state) phase diagram display non-trivial MCE with a se- 
quence of cooling and heating. 

Very recently, exact as well as numerical descriptions of 
the MCE in various one-dimensional interacting spin sys- 
tems have been attracted much attention 8-20 . Some two- 
dimensional systems have been also investigated, mainly, 
in the context of effect of frustration on the MCE 21,22 . 
The main features of MCE which have been revealed dur- 
ing the investigation of various models are: (i) essential 
enhancement of MCE in the vicinity of QCP, (ii) en- 
hancement of MCE by frustration, (iii) appearance of 
the sequence of cooling and heating stages during adi- 
abatic (de) magnetization for the systems demonstrating 
several magnetically ordered ground states, and (iv) po- 
tential application of MCE data for the investigation of 
critical properties of the system at hand. 

In this paper we continue the investigation of the MCE 
in one-dimensional quantum spin systems admitting the 
exact solution in the form of free spinless fermions via 
the Jordan- Wigner transformation (see, for instance, 
Ref. 23). Though, the cases of spin-1/2 XX (isotropic) 
and XY (anisotropic) models have been considered in 
Ref. 8, there is a series of spin chains with multiple spin 
interactions introduced by Suzuki in 70s 24,25 which can 
be solved by the standard Jordan- Wigner transforma- 
tion. We consider the simplest model of the Suzuki series, 
the spin-1/2 XX chain with three-spin interactions both 
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of XZX + YZY and XZY - YZX type 26 " 33 . As has 
been shown in previous investigations, inclusion of the 
three-spin interactions leads to appearance of new phases 
in the ground state and, thus, to more reach behavior 
in the vicinity of quantum phase transitions. We study 
two different types of three-spin interactions, namely, the 
three-spin interaction of the XZX + YZY type 27,28 and 
of the XZY - YZX type 26 ' 29 . Although these types of 
three-spin interactions are connected to each other by a 
unitary transformation 31 , in the pure form they repre- 
sent the systems with different symmetries and have dif- 
ferent ground-state phase diagrams. In particular, in the 
former case the ground-state phase diagram contains a 
point, where three different ground states merge [quan- 
tum triple point (QTP)]. Appearance of additional pa- 
rameters in the system, the three-spin coupling constants 
in our case, makes possible manipulation of the phys- 
ical features of MCE, namely, the position of the QCP 
and the values of the maximal and minimal temperatures 
during the adiabatic (de)magnetization. The knowledge 
about manipulation of the MCE physical parameters can 
be very useful for the future quest for the novel mag- 
netic materials and their applications in various aspects. 
On the other hand, the appearance of the points where 
several magnetically ordered ground states merge on the 
ground-state phase diagram, caused by the inclusion of 
additional three-spin interactions into the Hamiltonian, 
can lead to essential enhancement of MCE due to large 
entropy accumulation in such points. Finally, in real- 
life materials randomness is always present. It can be 
modeled assuming that on-site fields or intersitc interac- 
tions acquire random values. The considered quantum 
spin chains admit an exact analytical solution for ther- 
modynamics in the case when the transverse magnetic 
field is a random variable with the Lorentzian probabil- 
ity distribution 34 . As a result, with such a model it is 
possible to discuss the MCE in the presence of random- 
ness. 

The paper is organized as follows. At first, we present 
a general consideration based on the Jordan- Wigncr 
fermionization (Sec. II). Next, we consider separately 
the case of the of XZX + YZY interaction and the case 
of the XZY - YZX interaction (Sees. Ill and IV). Af- 
ter that, we consider a random-field spin- 1/2 XX chain 
with three-spin interactions (Sec. V). We discuss the 
MCE in all these cases. Finally, we draw some conclu- 
sions (Sec. VI). 



II. JORDAN- WIGNER FERMIONIZATION 
AND THERMODYNAMIC QUANTITIES 

Let us define the model under consideration. We con- 
sider N — > oo spins 1/2 placed on a simple chain. The 



Hamiltonian of the model looks as follows: 

JV 

« = zZ[-hs z n + J(s x n s x n+1 +slsl +1 ) 

+ K(s n s n+1 s n+2 + s^ l s n+1 s^ l+2 ) 

+ E{s x l sl l+1 s v n+2 — s^s^+isJi+2)] j (2-1) 

where h is the external (transverse) magnetic field, J is 
the isotropic XY (i.e., XX) exchange interaction con- 
stant (in what follows we will set J = 1 to fix the units), 
and K and E are the constants of the two types of three- 
spin exchange interactions. We imply periodic boundary 
conditions in Eq. (2.1) for convenience. 

The Hamiltonian (2.1) can be brought to the diago- 
nal Fermi-form after applying at first the Jordan- Wigncr 
transformation to spinless fermions, 

s+ — ? x +7's y — P ir* s~ — <i x — is y — P ,r 

c n Pn — l s n J c 'i = E'n — lS n , 

m m 

P m = H(l-2c f jCj ) = l[(-2s*), (2.2) 
and performing further the Fourier transformation, 

J _ p ikn J -ikn 

JV iV 
VJV n=l VJV n=l 

k = 2Trm/N, m = -N/2, N/2 - 1 (we assume that 
N is even without loss of generality). As a result, 

n = E**(4*-i). 

K E 
e k = -h + Jcosk cos(2fc) sin(2fc). (2.4) 



Using Eq. (2.4) we can easily calculate the partition 
function for the spin model (2.1) 

Z(T,h,N) = Tre-«/ T = II 2ch || ( 2 ' 5 ) 

k 

(we set &b = !)• Various thermodynamic quantities, such 
as the Hclmholtz free energy, the entropy, and the specific 
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heat (per site) immediately follow from Eq. (2.5): 
T\nZ(T,h,N) 



f(T,h) 



s(T,h) 



lim 

N->oo 



N 



1 

2^ 



dk ( — + T In n k 



c(T,h) = T 



df(T, h) 
dT 

1 

2- 

ds(T, h) 



T 



dT 



1 



2ttT 2 



dke k n k (l - n k ) 



(2.6) 



here n k = l/(e £fc / T + 1) are the occupation numbers of 
spinless fermions. Furthermore, we get 
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It may be useful to rewrite the formulas for thermo- 
dynamic quantities (2.6), (2.7) in terms of the density of 
states 

We have 
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(2.9) 



Formulas (2.9) are extremely useful for consideration of 
the random spin-1/2 XX chains within the Green's func- 



tions approach 34,35 , since that method permits to cal- 
culate the random-averaged density of states (2.8), see 
Rcf. 34 and Sec. V. 

Although the presented above formulas give a com- 
prehensive description of the quantum spin system (2.1) 
(and, in particular, the MCE), the thermodynamic be- 
havior is somehow hidden behind one-fold integrals in 
Eqs. (2.6), (2.7). More explicit dependencies of ther- 
modynamic quantities on temperature and field can 
be derived, e.g., in the low-tempcrature limit. Let 
us briefly discuss what happens with Eqs. (2.6), (2.7) 
when T — > 0, see also Sec. III. We note that n k (l — 
n k ) = l/{4ch 2 [e fe /(2T)]} and therefore as T -> 
only a small region where e k ~ is relevant in the 
integrals yielding c(T,h) or dm(T,h) / dT in the low- 
tempcrature limit. Clearly, if the energy spectrum of 
spinless fermions is gapped we immediately get that 
c(T, h) and dm(T, h)/dT vanishes as T -> 0. We turn to 
the case of a gaplcss energy spectrum of spinless fermions. 
Assume that we have e k = ef (k — h) z / ' z\ + . . . around 
k{ satisfying e ki = 0. Then we immediately find that 
c{T,h) oc T 1 ^' Note also that s(T,h) cx T 1 ^ and in 
a "flat-band-like" limit z — > oo the entropy becomes in- 
dependent on temperature (for a discussion of true flat- 
band spin systems see Refs. 8,13,14,36). While estimat- 
ing dm(T,h)/dT for odd z (e.g., for z = 1) we have to 
take higher-order terms in the expansion of e k around 
k = fcj. For even z [z = 2 for the QCP and z = 4 for the 
QTP, see Eq. (3.1)] we get dm(T,h)/dT cx T 1 '^ 1 and 
therefore T h cx T _1 . 

Alternatively the critical behavior can be derived us- 
ing formulas (2.9). The factor l/ch 2 [cj/(2T)] in the in- 
tegrands for c(T,h) and dm(T,h)/dT implies that in 
the limit T — > only a small region where uj « is 
relevant. Around the QCP p{uS) cx uj~ x I 2 and there- 
fore c(T,h) cx T 1 / 2 , dm(T,h)/dT cx T, whereas around 
the QTP p(u>) cx w~ 3 / 4 (see Ref. 31) and as a result 
c(T, h) cx 7 11 / 4 , dm(T, h)/dT cx T~ 3 / 4 . 

In the case of randomness considered in Sec. V, van 
Hove peculiarities in the density of states are smeared 
out, p(u}) has a finite nonzero value for any u> and, in 
particular, p(uj) = p(0) + . . . around u = 0. As a result, 
c(T, h) cx T for a sufficiently low temperature T — > 
[an estimate of dm(T, h)/dT requires higher-order terms 
in the expansion of p(to) around cj = 0]. The boundaries 
between different ground-state phases disappear since the 
quantum phase transition transforms into a crossover 34 . 



III. THREE-SPIN INTERACTIONS OF 
XZX + YZY TYPE 



Now we consider the chain with two-spin interactions 
and three-spin interactions of XZX + YZY type 27,28 . 
In this case we put in Eq. (2.1) E = and the energy 
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FIG. 1: The ground-state phase diagram in the K — h plane 
of the model (2.1) with J = 1 and E — 0. The dark-gray 
regions correspond to the spin-liquid II phase, the light-gray 
region corresponds to the spin-liquid I phase, and the white re- 
gions correspond to the ferromagnetic phase. The lines h*(K) 
which separate different regions correspond to quantum phase 
transitions between different ground-state phases. 



spectrum of spinless fermions in Eq. (2.4) reads: 

K 

Ek = —h + J cos k — — cos(2fc). 



(3-1) 



The third term in (3.1) may lead to a new ground-state 
phase, the so-called spin- liquid II phase 28 . The phase di- 
agram in the ground state is shown in Fig. 1. While \K\ 
is less than 1/2 ( J = 1) there are only two ground-state 
phases: the spin-liquid I phase and the ferromagnetic 
phase. However, when \K\ > 1/2 there is in addition one 
more ground-state phase, the spin-liquid II phase. It is 
worthy noting that there are two special points (K = 1/2, 
h = 3/4 and K = —1/2, h = —3/4) on the ground-state 
phase diagram at which all three ground-state phases 
meet (QTPs). For further details see Refs. 28,31. 

The entropy of the spin system is a function of the 
temperature T, the magnetic field h, and of the parame- 
ter K [see Eq. (2.6) in which we have to use given in 
Eq. (3.1) (J = l)]. 

Now we turn to a discussion of the MCE in its classical 
interpretation as an adiabatic change of the temperature 
of the considered model under field variation. Grayscale 
plots of the temperature at constant entropy s(T,h) = 
0.05 in the K — h plane arc shown in the upper panel of 
Fig. 2. The lowest values of the temperature are around 
the QTP K = 1/2, h = 3/4. Two lower panels of Fig. 2 
supplement the upper one. In two lower panels of Fig. 2 
we show by thin broken lines the dependencies T(h) at 
fixed values of s = 0.05, 0.10, . . . , 0.60 for spin model (2.1) 
with J = 1, K = 0.5, 1.5, £7 = 0. 

According to Eq. (1.1), to discuss the MCE we may 
analyze alternatively an isothermal change of the entropy 
under field variation. In the upper panel of Fig. 3 we 
show grayscale plots of the entropy in the K — h plane for 
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FIG. 2: Upper panel: Grayscale plots of the temperature as 
it follows from the condition s(T, h) — 0.05 in the K — h plane 
for model (2.1) with J = 1 and E = 0. Two lower panels: 
Isentropic dependence T vs h at s = 0.05, 0.10, . . . , 0.60 (from 
bottom to top in each panel) for model (2.1) with J = 1, 
K = 0.5, 1.5, E = 0. Thin broken lines correspond to the 
nonrandom model, thick solid lines correspond to the random- 
field model (5.1) with T = 0.1. 



constant temperature T = 0.05. Again two lower panels 
of Fig. 3 supplement the upper one. In two lower panels 
of Fig. 3 we show by thin broken lines the dependencies 
s(h) at T = 0.05, 0.10, . . . , 0.60 for spin model (2.1) with 
J = 1, K = 0.5, 1.5, £7 = 0. 

Comparing Fig. 2 and Fig. 3 with the ground-state 
phase diagram in Fig. 1 one can note that the MCE in the 
low-,s or low-T regimes perfectly reproduces the ground- 
state phase transition lines. 
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FIG. 3: Upper panel: Grayscale plots of the entropy s(T, h) 
in the K — h plane at T = 0.05 for model (2.1) with J = 1 and 
£7 = 0. Two lower panels: Isothermal dependence s vs h at 
T = 0.05, 0.10, . . . , 0.60 (from bottom to top in each panel) for 
model (2.1) with J = 1, K = 0.5, 1.5, E = 0. Thin broken 
lines correspond to the nonrandom model, thick solid lines 
correspond to the random-field model (5.1) with T = 0.1. 



Furthermore, consider, e.g., K = 1/2 (thin broken lines 
in the middle panel in Fig. 2). Clearly if we decrease 
adiabatically h from 2 to 3/4 the temperature noticeably 
falls down (e.g., approximately from 0.3388 to 0.0002 at 
s = 0.05 or from 0.4331 to 0.0025 at s = 0.10). We turn 
to the results shown by thin broken lines in the middle 
panel in Fig. 3. If we decrease isothermally h from 2 
to 3/4, the entropy of spin system noticeably increases 
(e.g., approximately from 0.00000 to 0.2195 at T = 0.05 
or from 0.00001 to 0.2670 at T = 0.10) meaning that 




-2.5 -2 -1.5 -1 -0.5 0.5 1 1.5 2 2.5 

h 

FIG. 4: T h vs h for model (2.1) with J = 1, K = (dotted), 
K = 0.5 (dashed), K = 1.5 (solid), and E = at T = 0.05. 
Thin lines correspond to the nonrandom model, thick lines 
correspond to the random-field model (5.1) with T = 0.1. 



the spin system absorbs from the thermostat the heat 
T [s(h = 3/4) - s(h = 2)] (that is, ps 0.0110 at T = 0.05 
or ps 0.0267 at T = 0.10) per site. 

In Fig. 4 we plot by thin lines the dependence Yh{h) 
at T = 0.05 for a few values of K = (dotted), K = 0.5 
(dashed), and K = 1.5 (solid). From the first glance it 
becomes clear that one peak (around h ss 0.75) is higher 
than others. It might be interesting to compare the 
height of maxima in the dependence Th(h) for T = 0.05, 
in particular the height of the high-field ones which cor- 
respond to cooling while h decreases starting from high 
fields (see thin lines in Fig. 4 for h = 0.75... 2). For 
K = 0, 0.5, 1.5 we have T h (h) ~ 11.00, 14.22, 10.80 
at h ps 1.02, 0.77, 0.93, respectively. (We recall 
here that the quantum phase transition occurs at h* = 
1, 3/4, 11/12 for K = 0, 0.5, 1.5, respectively, see 
Fig. 1 and Ref. 31.) Clearly, the value of the cooling rate 
around the QTP is about 130% of such a value around 
the QCP (e.g., for the studied earlier K = case 8 ). At 
lower temperatures the heights of maxima increase (e.g., 
T h (h) ss 55.56, 73.04, 55.41 at h ss 1.003, 0.754, 0.920 
for T = 0.01 and T h (h) ps 556.51, 739.48, 556.61 at 
h ps 1.0003, 0.7504, 0.9170 for T = 0.001) but the 
relation between the heights changes only very slightly 
and approaches roughly 3:4:3. Thus, one gets about 
33% larger change in temperature for the same adiabatic 
change of the field being performed around QTP in com- 
parison with being performed around QCP. 

This can be illustrated further while considering 
Eq. (2.7) in the limit T -> 0. We recall that accord- 
ing to Eq. (3.1) e k = —h + J — K/2 - (J/2 - K)k 2 + (J — 
8if)fc 4 /4! + . . . and therefore while approaching a high- 
field peculiar point h* from above, i.e., by decreasing h, 
h — > h* + (a ferromagnctic-to-spin- liquid transition), 
we have 

e k = -{h-h*)- 1 -e + ^ + ... (3.2) 
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for K = (QCP) with ft* = 1 and 

e k = -(ft - ft*) - -k 4 + —k 6 - —ft 8 
v ; 8 48 640 



(3.3) 



for K = 1/2 (QTP) with ft* = 3/4. Moreover, we may 
single out two regimes. In the first regime, we first take 
the limit T -> and then ft- ft* -> +0 [i.e., (ft-ft*)/T > 
1] . In the second regime, we first put ft — ft* — Q and then 
take the limit T -> [i.e., (ft-ft*)/T < 1]. [For the high- 
field peaks in Fig. 4 we have (ft - h*)/T 0.3 . . . O.4.] 
If T — > and (h-h*)/T > 1 wc can write n fc (l-n fc ) 
e — let I/t anc [ hence 



n fe (l-r, fc )oc e ^ 2 /^ 



for i-T = and 



n k {l - Hfc) oc e 



-fc 4 /(8T) 



(3.4) 



(3.5) 



for K = 1/2, see Eqs. (3.2) and (3.3). Furthermore, we 
can extend the limits of integration with respect to k 
in two relevant integrals in the formula for Th (2.7) to 
— oo and oo. Using Eqs. (3.2), (3.3), (3.4), and (3.5) we 
immediately find that Th — > l/(ft — ft*) as ft — s- ft* + in 
the limit T = for both cases, K = and K = 1/2. For 
small but finite T we obtain different results for K = 
and i-T = 1/2. Although all relevant integrals are doable 
with the help of the well-known formula for the gamma 
function T(z) (see, e.g., Ref. 37) 



I A -/V« r ( i_) ( 3. 6) 
a \a J 



9? A > 0, a, P > 0, it is simpler to find required results 
by using MAPLE codes. Namely, for the case K = wc 
have 



1 l + e/2-eT/8 
ft -ft* 1 + e- (ft-4)e 2 /4 



(3.7) 



with e = T/(h — h*), h* = 1, whereas for the case -K" = 
1/2 we have 



1 



1 + e/4- 0.119eVT + eT/32 



ft 



ft* 1 + e/2- 0.239eVT + (ft + 17/4)e 2 /16 



{3.8) 



with e = T/(h— ft*), ft* = 3/4. Approximate analytical 
formulas (3.7) and (3.8) yield 1 : 1.39 for ratio of the 
heights of peaks in the dependence T^(ft) around QCP 
and QTP at T = 0.05 that is in a reasonable agreement 
with exact numerical calculation according to Eq. (2.7). 

If we put at first ft = ft* and then assume T — > we 
can write 



integration with respect to k in the two relevant integrals 
in Eq. (2.7) can be extended to — oo and oo. After simple 
calculations using (3.6) we find that at QTP relates 
to T h at QCP as r(5/8)/r(9/8) to r(3/4)/r(5/4), i.e., as 
1.1267. . . :1. Interestingly, for arbitrary z we have Th — > 
{T[(z + l)/(2z)]/r[(2* + l)/(2z)]}/(2T) -* 0F/(2T) if 
z — > oo. Thus in such a case (z — > oo) Th is about 131% 
of T h at the QCP with z = 2. 

In summary, as can be seen from numerical calcula- 
tions of (2.7) reported in Fig. 4 as well as from an- 
alytical considerations in specific limits 38 , the efficiency 
of cooling while decreasing ft starting from the high-field 
limit is higher around the QTP than around the QCP. We 
note here that enhancement of the MCE in the frustrated 
J\ — J2 antifcrromagnctic Hciscnberg chain due to a can- 
cellation of the leading fc 2 -term in the one-particle energy 
spectrum at the point J2 = Ji/4 was discussed in Rcf. 8. 
For the case at hand (2.1), (3.1), the softening in the exci- 
tation spectrum occurs at the QTP which emerges owing 
to the three-spin interactions of XZX + YZY type. 



IV. THREE-SPIN INTERACTIONS OF 
XZY - YZX TYPE 

Now we turn to the spin-1/2 XX chain with three- 
spin interactions of XZY — YZX type 26,29 . In this case 
we put K = in Eq. (2.1) and the energy spectrum of 
spinless fcrmions in Eq. (2.4) reads: 



Ek = —ft + Jcosfc — — sin(2fc). 



(4.1) 



The ground-state phase diagram of the model is shown 
in Fig. 5. This type of three-spin interactions leads to 
the spin-liquid II phase too. But for the model (4.1), in 
contrast to the model (3.1), we do not have a QTP and 
all phase transition lines separate two different phases 
only, compare Fig. 5 and Fig. 1. For further details see 
Refs. 26,29. 

Figs. 6, 7, and 8 are similar to Figs. 2, 3, and 4. Again 
the MCE in the low-s and low-T regimes (Figs. 6, 7, 
8) indicates the ground-state phase transition lines seen 
in Fig. 5. Again we observe that cooling/heating is es- 
pecially efficient around QCPs. Since the model with 
three-spin interactions of XZY — YZX type does not 
have a QTP, in all cases Th behaves like it should around 
QCP, see thin lines in Fig. 8. For example, the height of 
the high-field peaks in the dependence r^(ft) at T = 0.05 
for various E is the same (see thin lines for ft — 1 ... 2 in 
Fig. 8). 



n k (l - n k ) 



2 + + e -e k /T 

I _ 1 -a 2 fc 2 V(4T 2 ) 



ik z ) 2 /T 2 



(3.9) 



with Ek = —ak z < 0, z = 2 for the QCP and z = 4 for 
the QTP, see Eqs. (3.2) and (3.3). Again the limits of 



RANDOM (LORENTZIAN) TRANSVERSE 
MAGNETIC FIELD 



In this section we use recent results on thermodynam- 
ics of the spin-1/2 XX chain with three-spin interactions 
in a random transverse field 34 to discuss the influence 




E 



FIG. 5: The ground-state phase diagram in the E — h plane 
of the model (2.1) with J = 1 and K — 0. The dark-gray 
regions correspond to the spin-liquid II phase, the light-gray 
region corresponds to the spin-liquid I phase, and the white re- 
gions correspond to the ferromagnetic phase. The lines h*(E) 
which separate different regions correspond to quantum phase 
transitions between different ground-state phases. 



of randomness on the MCE. To be specific, we consider 
the Hamiltonian (2.1) and make the change h — > h n , 
where h n is the random transverse magnetic field with 
the Lorcntzian probability distribution 



i 



r 



7T (h n - h) 2 + r 2 ' 



(5.1) 



Now h is the mean value of h n and the parameter T 
controls the strength of the Lorentzian disorder. The 
nonrandom case can be reproduced if T is sent to 0. All 
(random-averaged) thermodynamic quantities of the ran- 
dom quantum spin system (2.1), (5.1) can be expressed 
through the (random-averaged) density of states 34 . For 
instance, in the formula for the entropy in Eq. (2.9) we 
now have to use the density of states 

1 



p{w) 



GJM = - 



4 
A" 



(Zl - Z 2 )(z 1 - Z 3 )(Zl - Zi) 
~~2 



(Z2 - zx)(z 2 - z 3 )(z 2 - Zi) 



(5.2) 



and \z\\ < \z2 | < \za\ < \z^\ are the solutions of a certain 
quartic equation, sec Rcf. 34. In the case E = the 
relevant quartic equation, 



2J 



4 2 J 

z 3 + — (uj + h±iT)z 2 - — 2 + 1 = 0, (5.3) 
K K 



can be reduced to the quadratic one that further simpli- 
fies calculations. In the case K = the relevant quartic 
equation, 



2iJ 



4i 

z 3 + — (uj + h±iT)z 2 
E 



2iJ 



1 = 0, (5.4) 



£=1 




FIG. 6: Upper panel: Grayscale plots of the temperature as 
it follows from the condition s(T, h) — 0.05 in the E — h plane 
for model (2.1) with J = 1 and K — 0. Two lower panels: 
Isentropic dependence T vs h at s = 0.05, 0.10, . . . , 0.60 (from 
bottom to top in each panel) for model (2.1) with J = 1, K = 
0, E = 1, 2. Thin broken lines correspond to the nonrandom 
model, thick solid lines correspond to the random-field model 
(5.1) with T = 0.1. 



can be easily solved numerically. 

Our findings for the random- field models (2.1), (5.1) 
with r = 0.1 are shown by thick lines in the two lower 
panels of Figs. 2, 3, 6, 7 and in Figs. 4, 8. From the analy- 
sis of nonrandom models in Sees. Ill and IV we know that 
essential enhancement of MCE occurs around QCPs and 
QTPs. From the reported results for the random-field 
chains (compare, e.g., thin and thick lines in Figs. 4 and 
8) we conclude that just around these special points the 
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0.5 1 1.5 2 

h 

FIG. 7: Upper panel: Grayscale plots of the entropy s(T, h) 
in the E — h plane at T = 0.05 for model (2.1) with J = 1 and 
K — 0. Two lower panels: Isothermal dependencies s vs ft at 
T = 0.05, 0.10, ... , 0.60 (from bottom to top in each panel) 
for model (2.1) with J = 1, K = 0, E = 1, 2. Thin broken 
lines correspond to the nonrandom model, thick solid lines 
correspond to the random-field model (5.1) with T = 0.1. 

MCE is extremely sensitive to randomness. As can be 
seen in Figs. 4 and 8, even small randomness leads to a 
rounding of sparks in the dependence Th(h) and notice- 
able diminishing of maximal values of Th- 



VI. CONCLUSIONS 

In this work we have studied the MCE for the spin- 
1/2 XX chain with the three-spin interactions of the 




-2.5 -2 -1.5 -1 -0.5 0.5 1 1.5 2 2.5 

h 



FIG. 8: T h vs h for model (2.1) with J = 1, K = 0, E = 
(dotted), E = 1 (dashed), and E = 2 (solid) at T = 0.05. 
Thin lines correspond to the nonrandom model, thick lines 
correspond to the random-field model (5.1) with T — 0.1. 



XZX + YZY and XZY - YZX types. The consid- 
ered models have more parameters (in addition to the 
magnetic field we can also vary the strength of three- 
spin interactions), they contain several lines of QCPs 
and QTPs, and manipulation with MCE becomes pos- 
sible. We have found that the quantum phase transition 
lines clearly manifest themselves for the MCE in the low- 
s or low-T regimes. The ground-state phase diagrams 
can be perfectly reproduced by measuring of Th in the 
limit T -> 0. The vicinity of QCPs or QTPs is very ef- 
fective for cooling since low temperatures are achieved 
by only small decrease of the field. Particularly strong 
variation of temperature (at s = const) or of entropy (at 
T = const) with varying the magnetic field occurs in the 
vicinity of a QTP. 

We have discussed the MCE in a random quantum spin 
chain. We have found that even small randomness can 
noticeably diminish an enhanced MCE in proximity to a 
QCP/QTP. 

The considered models, thanks to their simplicity, 
have enabled the rigorous analysis of the thermodynamic 
quantities of interest. Although we do not know any par- 
ticular compound which can be described by the studied 
models, our results might have a general merit being use- 
ful for understanding the effects of proximity to QTP and 
randomness on the MCE. On the other hand, with fur- 
ther progress in material sciences and synthesis of new 
magnetic chain compounds, the lack of experimental data 
and comparison between theory and experiment may be 
resolved in future. 
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(2) / 1(1) 



r h -> -- 



E.VK 



(i) 



Alternatively, we may begin with P^ given in Eq. (2.9). As 
was mentioned in Sec. II, l/ch 2 [oj/(2T)] tends to 4TS(u>) 
in the limit T — s> and hence only a small region around 
ui — is relevant. Expanding p(w) around co — 0, p(uj) = 
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p(0) +p {1) (0)oj + 



we arrive at 
dp(u>) / dui\ u =o 

P(w)| u =0 



These two results for r\ are equivalent. 



